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ABSTRACT 

For the past decade, ionized outflows of a few 100 km s~^ from nearby Seyfert 
galaxies have been studied in great detail using high resolution X-ray absorption 
spectra. A recurring feature of these outflows is their broad ionization distri- 
bution including essentially ions (e.g., of Fe) from neutral to fully ionized. The 
absorption measure distribution {AMD) is defined as the distribution of column 
density with ionization parameter \dNH/d(\og^)\. AMDs of Seyfert outflows 
can span up to five orders of magnitude in C,- We present the AMD of five 
outflows and show that they are all rather fiat, perhaps slightly rising towards 
high ionization. More quantitatively, a power-law fit for log AMD oc (log^)" 
yields slopes of < a < 0.4. These slopes tightly constrain the density profiles 
of the wind, which until now could be addressed only by theory. If the wind is 
distributed on large scales, the measured slopes imply a generic density radial 
profile of n oc r~" with 1 < a < 1.3. This scaling rules out a mass conserving 
radial fiow of n oc r~^, or a constant density absorber, but is consistent with 
a non-spherical MHD outfiow model in which n oc r^^ along any given line of 
sight. On the other hand, if ionization variations are a result of local {6r) density 
gradients, e.g. as in the turbulent interstellar medium (ISM), the AMD slopes 
imply density scaling of n oc ^r""" with 0.7 < a < 1.0, which is quite different 
from the scaling of approximately n oc 6r^'^ found in the Milky Way ISM and 
typical of incompressible turbulence. 

Subject headings: X-rays: galaxies — galaxies: active — galaxies: individ- 
ual (IRAS 13349+2438, NGC3783, NGC 7469, NGC5548, MCG -6-30-15) — tech- 
niques: spectroscopic — galaxies: ISM 
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1. INTRODUCTION 



The X-ray spectra of many active galactic nuclei (AGNs) viewed directly toward the 
central source (e.g., Seyfert 1 galaxies) show the continuum flux absorbed by numerous 
absorption lines produced by ionized gas. The role of these absorbers for the AGN and for the 
host galaxy are still being debated. In order to advance our understanding of these outflows 
and their significance to the AGN, it is important to develop model-independent techniques 
for quantifying the physical parameters of the wind. So far, only the wind velocity measured 
from Doppler shifts and the total column density Nh deduced from the equivalent widths 
of absorption lines are known with high confidence. The broad range of ionization (and 
temperature) in the wind is also unambiguous. Conversely, there is httle robust constraints 
on the location, geometry, and density of the wind, which leaves its mass outflow rate, 
momentum, and energy fairly uncertain. 

The X-ray band can potentially provide the full physical picture of AGN outflows, as 
it comprises detectable absorption lines from the full range of charge states: From neutral 
up to hydrogen-like. The detection of all of the absorbing charge states of a given element 
enables an accurate measurement of the total column density. For this purpose, Fe ions are 
particularly useful as they form over several orders of magnitude of ionization parameter ^ 



'V. 



where L is the ionizing luminosity between 1 and 1000 Rydberg, n denotes the hydrogen 
number density and r is the distance of the absorber from the ionizi ng source. AGN outflows 



can be distributed over the full range of ionizatio n — 1 < log,^ < 4 (ISteenbrugge et al.l 12003 



2005 



tion. 



Holczer et al 



Costantini et al.l 120071 : iHolczer et al.l 120071 ). In order to best quantify this distribu- 



(120071 ) defined an absorption measure distribution [AMD), which is the 



distribution of the hydrogen column density Nh along the line of sight as a function of log ^ 



AMD = \dNH/d{\og^)\ (2) 

Reciprocally, the total Nh can be expressed as an integral over the AMD. The AMD 
is the absorption analog of the Emission Measure Distribution (EMD) widely used in the 
analysis of emission-line spectra. It provides a more complete representation of the ionization 
distribution than the more commonly used models of several ionization components each 
with a fixed ^; the overall success of the latter in producing good spectral fi ts to the data 



notwithstanding. The AMD reconstruction method is outlined in detail in iHolczer et al. 



(120071 ). Here, we only note that for each target, the observed continuum spectrum is used to 
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obtai n the ionization balance as a function of from the XSTAR code (iKallman &: Bautista 



200ll ). Subsequently, the AMD is obtained by a fitting process that aims to best reproduce 



all of the individually measured ionic column densities. Being the solution of an inversion 
problem, the reconstructed AMD is inevitably degenerate. Nonetheless, a step function is 
the preferred AMD form, as it allows rigorous local error calc ulations that t a ke int o account 
the degeneracy between different ^ bins. For more details see iHolczer et al.l (120071 ). 

In this paper, we present the available AMD distributions, which were derived for five 
sources with high signal-to-noise-ratio X-ray grating spectra. The objective is not to focus 
on any individual target, but to seek trends in the best available (yet admittedly small 
and incomplete) sample of AGN outfiows. The idea is to link the general AMD behavior 
with the physical properties of the outfiows, in particular the density profile, in order to 
obtain an insight into the origin and astrophysical significance of the winds. These physical 
implications are described in ^and discussed in §1] 



2. Targets and AMDs 



We include all the Seyfert outfiow targets with sufficiently good grating spectra for 
which we were able to reconstruct a non-degenerate AMD. The AMDs o f IRAS 13349+2438 
and of NGC 3783 were presented in great detail in iHolczer et al.l (120071 ) and are based on 
deep exposures with the HETGS spectrometer on board Chandra. The AMD of MCG-6- 
30-15 is from a recent analysis of the HETGS spectrum (Holczer et al., in preparation) and 
corresponds to the slow AGN component in that complex spe ctrum. The AMD o f NGC 5548 
was derived from the ionic column densities published by ISteenbrugge et al.l (120051 ) and 
based on spectra from the HETGS and LETG instrumen ts on board Chandr a. The AMD 
of NGC 7469 is an improved version of that presented in iBlustin et al.l (120071 ) and is based 



on spectra acqui red with the RGS sp ectrometer on board XMM-Newton. The improvement 
i s tha t unlike in iBlustin et al.l ( 20071 ). here we followed the general method of IHolczer et al. 



( 120071 ) and allowed the AMD to be zero in certain bins, if required by the fit. We note that a 
th orough spectral analysis was carried out also for the Chandra LETG spectrum of Mrk 279 
by I Cost ant ini et al.l (120071 ). However, the low column density and marginal detection of Fe 
absorption in Mrk 279 is insufficient to produce a well constrained AMD. 

The five AMDs used in the present work are plotted in Fig. [H Note that all AMDs 
have a gap at intermediate ^ values, which is in fact a data point at zero that is not plotted. 
For four sources this gap fall s between log £ ^ 0.5 — 1.5 ergs~^cm, wh il e for NGC 7469 it is 



at slightly higher values. In IHolczer et al.l (120071 ) and in IBlustin et al.l (120071 ). these AMD 
discontinuities were interpreted as regions of thermal instability. A physically separated two- 
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component outflow is another possible interpretation, although the similar outflow velocities 
of the low and high ionization regions and the ubiquity of the minimum at approximately 
the same temperature regime argue against two separate components. The presence of the 
thermally unstable region and its exact position (in C.) dep end strongly on the corresp ondence 
between charge state and ^ that are derived from XSTAR ( iKallman fc Bautistall200ll ) and are 
uncertain to some degree. It is known, for instance, that some of the available die l ectronic 
recombination (DR) rates used in these computations are inadequate (jNetzerll2004j : iBadnell 



20061 ) . Another effect that could potentially change the unstable region by affecting the 



^ at which low charge states form in the outflow is photoionization of these ions by the 
EUV continuum that is poorly constrained. Finally, all ionization balance computations are 
carried out with a single ionizing spectrum, i.e., in the optically thin limit. This is justified 
by the observed X-ray spectra that typically show only shallow photoelectric edges, if any. 
Once the medium becomes optically thick and deeper regions are exposed to lower (and 
harder) ionizing luminosity, this approximation gradually loses its validity. 

In this work, we are interested more in the broad characteristics of the distribution. 
Assuming that the minimum is due to thermal instability justifies viewing the remaining 
stable regions of the AMD as different ionization components of the same outflow that are 
likely driven by the same physical mechanisms. We thus apply simple (least mean square) 
linear fits to each AMD according to log AMD{^) = alog^ + 6, or AMD{C,) = lO^C- 
The allegedly thermally unstable ionization regi ons are ignored in t he fit. A similar power 
law distribution of ionization was suggested by ISteenbrugge et al.l (120031 ). It can be seen 
in Fig. [T]that a power law representation of the AMD for some targets is far from perfect. 
Nevertheless, given the uncertainties on each data point, and for the mere purpose of seeking 
AMD trends, we deem it sufficient. Moreover, a power law AMD representation allows 
direct comparisons with self-similar wind models as discussed in the following sections. The 
best-fit power laws for the five outflows are shown in Fig. [T] and the fit parameters are listed 
in Table m Evidently, all five power laws have relatively similar slopes of a = 0.0 — 0.4. The 
implied density pr o file sl opes (a), introduced in the next section, are also hsted in Table [TJ 
Steenbrugge et al. J2005h fitted a power law to their distribution of individual ionic column 
densities in NGC 55 48, which is somewhat similar but not identical to the full AMD (see 
Holczer et al.l 120071 . for details). They obtained a slope of a = 0.2 ± 0.1 (denoted there a) 
with Fe ions alone, which is in excellent agreement with our best-fit of a = 0.14 ±0.08. When 
they use other elements, however, and need to rely on the assumption of solar abundances, 
they get a somewhat steeper slope of a = 0.4. 
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3. Physical Implications for Outflows 

In this section, we explore what the AMD shape, namely its slope, implies for the 
density profile of the outflow. It is convenient for this purpose to consider two separate 
boundary cases. One is a radial large-scale outflow in which the density n(r) varies as a 
function of the distance from the source r. In this case, ^ = L/{nr'^) varies along the flow 
due to the combined change in r and in n. The other possibility is a remote absorber at a 
distance more or less constant r = ro, but with strong density gradients over small distances 
5r <^ tq. In this case, the ionization varies exclusively due to the local density gradients. 
Obviously, a real outflow could be a combination of the two, namely a large scale outflow 
with strong local density gradients. However, for the current analysis, it is instructive to 
consider these two cases separately. 



3.1. Large Scale Outflow 

A large-scale outflow is deflned here as an outflow with a size comparable to the distance 
r from the ionizing source so that the ionization parameter (eq. [T]) is a function of r also 
through the parameterization of n(r). Let us assume a power-law density proflle with an 
index a > 



n{r) oc r ° (3) 

Such a parameterization of the density is convenient for comparison with self-similar wind 
solutions. The ionization parameter thus varies with radius as ^(r) oc r"~^, or reciprocally 

r{i) oc ^V("-2) and 

— oce -2 (4) 

The radial increment of hydrogen column density with the ionization parameter ^(r), 
can be written as 

dNn = n{r)dr = n{r)^d^ (X ^~'^d^ (5) 
Consequently, the AMD (eq. [2]) scales with ^ as 
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AMD 



dN, 



H 



dhg^ 



H 



OCT 



g-l 
' Q-2 



(6) 



where a = — (a — l)/(a — 2) is precisely the slope of log AMD vs. log^ fitted for in Fig. [H 
Finally, a and its uncertainty can be expresses as 



l + 2a , 
1 + a 



Aa 



[l + a) 



(7) 



where Aa stands for the standard error of the fitted power law slope a. The resulting a 
values for the five outflows are listed in Table [TJ Note that a scaling with ^ similar to that of 
eq. is obtained for the integrated Nh up to r, as Nnir) oc n{r')dr' oc r^~" oc . 
However, Nh{t) is a property that pertains to the entire flow up to radius r, while the 
AMD is a local property of the flow at each position (or as a function of ^), and is therefore 
preferred here. 

Note that for AMDs that increase with ^ (i.e., flow dominated by high charge states) as 
found in the present sample, a must obtain values between 1-2. The observed slopes, e.g., 
of a = 0, 0.25, or 0.4, imply density profiles that are not particularly steep with respective 
indices of a = 1.0, 1.2, and 1.3. A steep AMD over narrow ranges of ^, on the other hand, 
is difficult to rule out. Indeed, by inspecting Fig. [H it can be seen that some outflows (e.g., 
NGC3783, MCG -6-30-15, NGC5548) tentatively show steeper AMDs at low ionization (^ 
that flatten out at higher ^. By analyzing eq. ([6]), a few special cases can be easily examined. 

For a constant- density Uq absorber, a = 0, and the AMD would slowly decrease with 
ionization parameter as AMDino) oc i""^^"^. No slope with a < is observed in the five 
present outflows, which rules out a constant density outflow. Indeed, the constant density 
scenario is quite unphysical. 

Yet another possible wind scenario is a steady mass conserving radial flow similar to 
a stellar wind, in which the mass outflow rate, opening angle, and wind velocity are all 
constant. This results in a de n sity p rofile n oc r~^, or a = 2. Such a flow was assumed 



for instance by iKrolik &: Krisd (119951 ). In terms of the AMD slope a, it tends to infinity 
(eq. in]). In this case, the ionization parameter would of course be constant (^o) along the 
flow implying a sharply peaked AMD (around ^o)- This is in stark contradiction with the 
present sample that all feature broad AMDs. It is worth noting however, that some Seyferts 
do show distinctly faster (few 1000 km s~^) and more highly-ionized (log^ > 3) outflows 
that perhaps could be characterized by a single ionization parameter, a good example of 
which is the high-velocity components (4500 km s~^ & 1900 km s~^, respectively) found 
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in NGC 4051 flSteenbrugge et al.l l2009[ ) and in MCG -6-30-15 flSako et al.ll2003h . The low 
and intermediate charge states are clearly absent from these high-velocity components that 
do not have the same broad AMD as the present, slower winds. This along with the lack 
of a continuous velocity distribution (c.f., distinct slow and fast components) suggests that 
the slow and fast outflows are different types of winds, and may also differ in their position 
within the AGN system. 

All of the outflows observed to date for which an AMD could be reconstructed imply 
a fiat to modestly increasing AMD with ^. The column density measured for the low- 
ionization species is generally never higher than that of the high ionization plasma in terms 
of equivalent Nh- In order to obtain a strictly fiat AMD while retaining the assumption 
of a global power-law density dependence on radius (eq. [3]), the density scaling must obey 
n(r) oc 1/r (or a = 1). Allowing for a moderately increasing AMD with ^, as observed in 
Fig. [1] results in density scalings r~" with 1 < a < 1.3. See Tabled] for details. 



3.2. Small Scale Outflow 

Perhaps an equally interesting case is a distant absorber (at tq) in which the ionization 
parameter distribution ^ ^ L/{nrQ) is driven by density gradients on much smaller distances 
{6r <^ ro) . The broad AMD then requires the density to vary by several orders of magnitude 
to account for the different orders of magnitude of observed ^. Such a multi-phase scenario 
could be a natural consequence, e.g., of turbulence as observed in the interstellar medium 
(ISM) gas in our ga l axy. I ndeed, such a description for Seyfert outflows was proposed by 
Chelouche fc Netzerl (200511 for NGC 3783 and further systematically explored for similar 



outflows in 



Chelouchd fl2008f ). 



We now assume the density varies locally on these small length scales 5r 



n{5r) oc 5r 



which merely reflects typical density variations with distance that are assumed to be isotropic. 
In this case, the dependence of the AMD on a is slightly different than in Eqs. (jl]-[6]), as 

^ oc n"-*^ oc (5r", subsequently ^^^^ oc , dNu = n{6r)^^^^^d^ oc ^~^^d^, and the AMD 

scales with ^ as 



AMD{r = ro) 



dN 



H 



dXogi 



dN, 



H 



di 



(9) 
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where again a is the AMD slope and thus here a(r = tq) = —{a — 1)/q;, and inversely 



a 



1 Aa 
± 



'1 



(10) 



which can be compared with the large-scale flow relation given in eq. ([7]). Here again, Aa is 
the standard error of the fltted power law slope a. The resulting values for the flve outflows 
considered here are listed in Table [TJ The functional forms of the density scaling a in eqs. 
([7]) and ffTOl) are generally different. However, in the relevant regime of observed slopes a ^ 
(and only there), the value of a implied by both expressions is similar and close to unity. We 
obtain 1 < a < 1.3 for the large scale flows and 0.7 < a < 1 for the small scale (constant 
ro) flows. Far from a = 0, the two would be signiflcantly different. 



It has been shown by iGoncalves et al.l (120061 ) that strong density gradients can occur 

jalanc e as thermal instabihties (suggested 
' 2007h result in strong temperature jumps 



naturally in an absorbing slab that is in pressure 



also by the discontinuous AMDs in iHolczer et al. 



that in turn cause strong density gradients. IChelouche &: Netzerl ( 120051 ) propose a different 
model for NGC 3783 of thermal launching of a wind with assumed strong density gradients on 
small length scales and invoke a density scaling similar to eq. ([8]). They obtained a density 
profile with a = 0.8 (deno ted there as —jS), which is in very good agreement with the present 
results. Chelouche (l2008h recen t ly app lied the same model to a sample of outflows. However, 
the values quoted in lChelouchd (120081 ) are a ~ 1.3 ± 0.6, systematically higher than what we 
flnd here, although still consistent given their large errors. In any event, the present method 
provides much tighter constraints on a (Tabled]). 



4. Discussion 

The AMD slope deflned in Eqs. ([6]) and (Q provides a novel way to obtain a self 
similar density proflle of the outflow. The slopes available for a handful of outflows are all 
quite similar and lie in the range of < a < 0.4. The fact that a is not much larger than 
zero reflects the broad and relatively flat ionization distribution of Seyfert outflows. Along 
with the moderate velocities of a few 100 km s~^, these are the only two model- independent 
recurring features of such flows. In ^ we quite generally differentiate between two possible 
physical flows that could produce this AMD proflle. One is a large scale AGN outflow that 
would need to have a density proflle of n oc r^" with power-law indices of 1 < a < 1.3 in 
order to reproduce the observed AMD. Importantly, a smooth radial mass conserving flow 
of n oc r~^, for example, that would correspond to a single or a sharply peaked AMD are 
conclusively ruled out, and so is a const ant- density absorber (a = —0.5). Another possibility 
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is a remote absorber with strong density gradients on local length scales 5r in which case 
n oc 5r~" and 0.7 < a < 1. 

If the outflow originates in the accretion disk, as believed by many authors, it is of the 
large scale type. The derived density profiles can thus be compare d directly with an entire 
family of self-similar solutions for hydromagnetic (MHD) disk winds (IContopoulos fc Lovelace 



Blandford k. Payii^ 



19941) . To obtain an MHD wind with a slope of a = 1.5, for example, as in 
( jl982 ) would require a global AMD slope of a = 1.0 that is clearly excluded by the data. 
Such a high index is not seen i n any of the Seyfert outflows (Fig. [H Table [T]). Models from 
Contopoulos &: Lovelacel (1994), on the other hand, with 1 < a < 1.3 are consistent with 
the observed AMDs (the present a corresponds in their n otation to 3 — 2a: ). In f act, the 
specific self-similar solution with a = 1 was preferred by iKonigl &: Kartjd (119941 ) for its 
minimal magnetic energy. It is interesting to note that a totally different approach using 
the temporal behavior of accreting sources suggests indices of a = 1 — 1.5 for the Comp- 
ton scattering rnedium that produces the continuum X-ray spectrum (IKazanas et al.l 119971 : 
Hua fc KazanasI 1997). Consequently, a detailed MHD wind model has been developed by 
Fukumura et al. (in preparation) with an n oc profile that by definition produces a fiat 
AMD. Note that density profiles th at significantly deviat e from n oc have been reported 
out to the narrow line region (NLR). iBennert et al.l (120061 ) measured values of a = 1.46 ±0.2 
and a = 1.14 ± 0.1, respectively, for a sample of Seyfert 1 and Seyfert 2 NLRs. 

For local density gradients, it is instructive to think of pressure-balanced clouds. This 
type of models uses an absorbing slab to represent the AGN outflow, and with careful numeri- 
cal treatment of the radiative transfer problem through the sla b can compute its temperature, 
density, and ionization distributions (IGoncalves et al.ll2006l ). The astrophysical context of 
su ch "slabs" and h o w the y are launched remains to be sorted out, but models such as those 
of IGoncalves et al.l (120061 ) already stand out by virtue of their continuous and broad AMD, 
which is observed for Seyfert outflows. In fact, with some fine tuning, these models can 
reproduce the majority of the measured ionic column densities, as well as the observed un- 
stable region of the AMD. On the other hand, they also predict other unstable regions that 
are not necessarily observed as such. The numerical nature of these radiative transfer codes 
does not straightforwardly provide for an analytical power law density profile to be directly 
compared with eqs. ([HD and (fTUl) . 

An alternative natural source of density gradients on many scales is turbulence, such as 
that found in the interstellar medium (ISM) of the Milky Way. For a recent re view of ISM 
observations and the relevant physical processes, see lElmegreen fc Scald (120041 ). A general 
feature of the ISM is its power spectrum 
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P{k) oc k- 



observed to be appro ximately a power law o ver five decades in wavenumber k = l/6r from 



10~ — 10~ cm~ ( lArmstrong et al.lll995l ). The standard definitions are used here for 
P{k) = n{k)n{k)* where fi{k) is the Fourier transform of an observable quantity (e.g., 
density) n{6r), and n{k)* is its complex conjugate. Typically, 3D ISM power spectra are 
consistent with, or slightly steeper than the Kolmogorov spectrum for non-magnetic, incom- 
pressible turbulence in which p = 11/3, as a result of energy ca scading from large to small 
scales with no preferred intermediate scale (IKolmogorovl 1 1 94 1 ) . So mewhat steeper power 
spectra {p = 4) can be expected from random shocks ( SaffmaiJ 1971 ). While in a magnetic, 
compr essible turbulent medium, d ensity fluctuations can scale differently on different length 
scales ( iLithwick fc Goldreichll200ll ). the prevalent spectrum over most scales is a Kolmogorov 
one. Hence, for the most part a value of p = 3.8 ± 0.2 is a reasonable approximation for the 
ISM. 

In order to obtain a density scaling akin to eq. ([H]) for turbulent media, the mean squared 
density variations (variance) can be estimated from the volume 3D integral in Fourier space 
of the power spectrum (Parseval's theorem) over a cell the size of l/5r 



{\5n^\)= / P{k)d^k oc 6rP- 



(12) 



The root mean square (rms) of the density fluctuations in a volume ~ 6r^ thus scales as 



p-3 



(13) 



If this rms of density fluctuations is understood as the typical scaling of density variations 
over Sr as in n{Sr) oc 6r~°' with a = (3 — p)/2, it can be directly compared with eq. ([8]). Note 
that since generally p > 3, the density fluctuations do not vanish at any scale within the 
medium (i.e., no asymptotic homogeneity). Thus, eq. (JISD, in principle, could be meaningful 
over many orders of magnitude. 

Applying the typical ISM value of p = 3.8 ± 0.2 yields a density scaling of oc 6r^'^ 
{a = — 0.4± 0.1). Recall that using the AMD slopes, we obtained markedly different density 
gradients with 0.7 < a < 1.0 (Tabled]). Another difference between the Seyfert outflow and 
typical ISM is the range in density. For a distant, localized (ro) Seyfert absorber, the range 
in density required for producing the observed ionization structure needs to span up to five 
orders of magnitude (— 1 < log^ < 4). It is not clear that the ISM spans such a broad density 
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range, and certainly not within a single medium (e.g., a molecular cloud). The observed five 
orders of magnitude in k in the ISM correspond to only (n oc two orders of magnitude 

in density. In order to obtain the high density gradients of 0.7 < a < 1.0, an extremely 
shallow power spectrum with p = 0.8 — 1.0 would be required. We are not aware of any 
other (ISM) medium with such a flat power spectrum. 

5. CONCLUSIONS 

We have presented the Absorption Measure Distribution (AMD) of five nearby AGN 
outflows. The AMD is the differential column density through the flow as a function of ion- 
ization parameter ^ that can be derived from the high resolution X-ray absorption spectrum 
of the outflow. Key to reconstructing the AMDs is the observation of almost all Fe charge 
states from neutral to H-like. Using a power law approximation, we flt for the AMD slopes 
and find that there is only a mild increase of column density with ^, namely the slopes are 
all in the range of to 0.4. It may be possible that at low ionization log^ < 1, the AMD 
is steeper and then fiattens out towards higher ionization (log^ > 2). However, the limited 
ionization resolution in the AMD (due to the large formation regions of the various charge 
states) does not allow for a more precise assessment. In any event, the AMD reconstruction 
at low ionization is currently hampered by the uncertain dielectronic recombination rates 
and by the unknown EUV continuum. We expect the former limitation to be remedied 
soon by new laboratory measurements and improved atomic calculations. We also expect 
more high-quality grating spectra of AGN outfiows to be released to the Chandra and XMM- 
Newton archives in the next few years. Future high throughput X-ray spectrometers, such 
as the SXS currently being built for Astro-H will further increase the sample. As the sample 
of well studied Seyfert outflows increases, the present method will continue to facilitate our 
growing understanding of the physics governing these winds. 

We have discussed the best-fit AMD slopes in the context of two qualitatively different 
scenarios, one is a continuous large-scale radial ffow and the other is a distant well localized 
absorber. If the outflow is of the former type, we can constrain its density proflle rather 
tightly to n oc r~" with 1 < a < 1.3. This already rules out the standardly assumed n oc 
proflle, or even a proflle of n oc r^^'^. A constant-density absorber (a = 0) is also clearly ruled 
out. An MHD wind, on the other hand, with a shallower density proflle of a > 1 as proposed 
by Fukumura ct al. (in preparation) is consistent with the observed outflows. Alternatively, 
slab models with proper radiative transfer can produce a broad ionization distribution, with 
intermediate thermally unstable regions, generally consistent with the observed AMDs. If 
the absorber is remote from the central source, and compact, then its ionization depends 
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solely on the local density gradients. We find that in this case, the density should vary 
as n oc (5r~° with 0.7 < a < 1. Such sharp gradients, and occurring over five orders of 
magnitude in density, are quite different from those typically found in the Milky Way ISM 
and interpreted as due to turbulence. 

I thank Tomer Holczer for providing the AMD results in simple data form and Alex 
Blustin for suggesting the title for the paper. 
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Table 1. Best-fit parameters from tlie linear fits in Fig. ^AMD(^) oc 



Target 


Slope (a) 


a ^ 


a'' 


NGC 3783 


0.29 ± 0.11 


1.22 ± 0.07 


0.78 ± 0.07 


IRAS 13349+2438 


0.02 ± 0.01 


1.02 ± 0.01 


0.98 ± 0.01 


MCG -6-30-15 


0.10 ± 0.11 


1.10 ± 0.09 


0.91 ± 0.09 


NGC 5548 


0.14 ± 0.08 


1.12 ± 0.06 


0.88 ± 0.06 


NGC 7469 


0.24 ± 0.13 


1.19 ± 0.08 


0.80 ± 0.08 



'^For a distance driven AMD, a = (1 + 2a)/(l + a), as in H3.1\ 
'Tor a density driven AMD, a = 1/(1 + a), as in 
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Fig. 1. — AMD distributions obtained for five targets. The best fit power-law is presented 
for each AMD. The inferred power-law slopes are relatively moderate, indicating a rather 
flat distribution of column density with ionization parameter. 



